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Abstract 

Monte Carlo simulations are used to show that the steady state of the d = 2, 
two-temperature, diffusive XY model displays a continuous phase transition 
from a homogeneous disordered phase to a phase with long-range order. The 
long-range order exists although both the dynamics and the interactions are 
local, thus indicating the failure of a naive extension of the Mermin- Wagner 
theorem to nonequilibrium steady states. It is argued that the ordering is due 
to effective dipole interactions generated by the nonequilibrium dynamics. 
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Fluctuations increase as the spatial dimension of a system is decreased and low- 
dimensional, equilibrium systems with finite-range interactions cannot maintain long-range 
order. This observation is known as the Landau-Peierls (LP) argument |IJ or, for systems 
with continuous symmetry, as the Mermin- Wagner (MW) theorem ||. Since phase transi- 
tions and emerging long-range order in nonequilibrium steady states have attracted a lot of 
interest recently ||||], it is natural to ask if similar arguments are valid for nonequilibrium 
systems. 

Clearly, the LP argument or the MW theorem must be generalized if we want to apply 
them to a nonequilibrium steady state. The reason is that while interactions and temperature 
define an equilibrium state, a nonequilibrium steady state is the result of the interplay of 
interactions and dynamics. Thus apart from restricting the range of interactions, we must 
also qualify the dynamics in a nonequilibrium process. To understand this point, consider 
a kinetic Ising model in which two dynamical processes compete with each other ||. The 
elementary processes are (i) spin flips generated by nearest-neighbor interactions coupled 
to a heat bath at temperature T and (ii) random spin exchanges generated by a T = oo 
heat bath. In this model, the interactions are short ranged but the dynamics is long- 
ranged because spins are exchanged at arbitrarily large distances. As a result of long-range 
dynamics, effective long-range interactions are generated and long-range order appears in this 
flip-and-exchange kinetic Ising model || even in d — 1. This result is in contradiction with 
the LP argument which states that long-range order is impossible for a d = 1 equilibrium 
Ising system with short-range interactions. Thus a simple extention of the LP argument 
and of the MW theorem to nonequilibrium systems without restricting the dynamics clearly 
doesn't work. 

One might conjecture that by restricting both the interactions and the elementary dy- 
namical processes to be local, the generation of effective long-range interactions would be 
prevented and an extention of the LP or MW arguments to nonequilibrium systems could 
be made. This conjecture, however, is not necessarily true since local interactions combined 
with anisotropic local dynamics have been shown to lead to effective long-range (dipole) 
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interactions in a two-temperature kinetic Ising model j7|. Furthermore, a spherical approxi- 
mation to the above model has shown || that effective long-range interactions are generated 
not only in the Ising case but also in the case of infinite-component order-parameter field. 
Thus one may expect that neither the LP argument nor the MW theorem can be extended 
to nonequilibrium steady states even if both the interactions and the dynamics are local. 
Of course, this expectation should be checked carefully since the details of the dynamics 
do play an important role in determining the ordering properties of nonequilibrium steady 
states, and the details, such as the presence or absence of topological defects, the type of 
defects determining the long-time relaxation of order, are indeed very different in the Ising 
model, in the spherical limit and in the case of an n-component order parameter field with 
1 < n < d. 

Our aim here is to confirm the failure of the naive extention of the MW argument by 
studying the existence of long-range order in the d = 2, two-temperature, diffusive kinetic 
XY model. According to the MW theorem, long-range order cannot exist in the d = 2 
equilibrium XY model with short range interactions. However, we will present the results of 
Monte Carlo simulations which show that a continuous transition from the homogeneous dis- 
ordered state to a state with long-range order exists in the nonequilibrium two-temperature 
model. We will also construct a field theoretic description of the transition that is consistent 
with the Monte Carlo results. The theory indicates that, just as in the two temperature 
Ising case HJ, the critical properties are in the universality class of an equilibrium system 
with long-range interactions, supporting our conjecture that the generation of effective long 
range interactions is the mechanism responsible for the breakdown of the MW theorem. 

The model we study consists of XY spins rhi, (i. e. of two dimensional vectors of unit 
magnitude) at sites i of a square lattice with periodic boundary conditions imposed. There 
are ferromagnetic interactions of strength J between nearest-neighbor spins, so the energy 
of a configuration is H = —JY^Uj)^ • ^j, where (ij) indicates nearest-neighbor pairs. 
The system evolves by Kawasaki exchanges |5] of nearest-neighbor spins, occuring with 
Metropolis rates flnj. Exchanges along one axis of the lattice, which we call the 'parallel' 
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direction satisfy detailed balance at an inverse temperature j3\\ = J/T\\, while exchanges in 
the other direction, which we call the 'perpendicular' direction satisfy detailed balance at an 
inverse temperature (3± = J/T±. When the temperatures of the two heat baths are equal, 
the dynamics satisfies detailed balance with the nearest neighbor XY Hamiltonian, and the 
model reduces to the equilibrium kinetic XY model with Kawasaki dynamics. However, 
when the two temperatures are not equal, an energy current flows from the high-T heat 
bath to the low-T heat bath, and detailed balance is not satisfied. The critical behavior of 
the system when the two temperatures are not equal is very interesting in the Ising version 
of the above model f7j. Its universality class coincides with that of another equilibrium 
model: The Ising ferromagnet with dipolar interactions . As we shall see below analogous 
situation develops in case of the XY model. 

For simplicity, in our simulations, we considered only the case /3u = 0, i.e. exchanges 
in the parallel direction were random. Furthermore, the distribution of the angles of the 
spins (which is conserved by the dynamics) was chosen to be uniform: for a system with 
iV = L±_ ■ L» lattice sites the angles of the N spins {?%} were 2itj /N with j = 1, 2, N. A 
difficulty in studying the phase transition in this system is that the spatial anisotropy requires 



an analysis using anisotropic finite-size scaling [IT|. That is, one must compare systems 
whose shapes scale such that L±/Ll +A is constant, where A is an anisotropy exponent. As 



was done for the two-temperature Ising system |12 |, and based on renormalization group 
results (discussed below) which indicate that A ps 1 in the present model, we choose to 
simulate systems with sizes L±_ x L\\ = 12 x 9, 16 x 16, 24 x 36, and 32 x 64, which are 
related by the naive scaling 4L x L 2 . 

The simulations revealed a homogeneous disordered state at small /3±, and an ordered 
state with long-range order at large /3±. A typical ordered configuration is shown in Fig. 1. 
Because the spin distribution is conserved by the dynamics, the ordering occurs as a phase 
separation resulting in a steady-state configuration with a spin wave in the perpendicular 
direction. To study the transition with the type of ordering shown in Fig. 1, we defined the 
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order parameter, as the following average of long- wavelength limits of structure factors 

* = -[C^tt/Ll, 0) + C 2 (2n/L ± , 0)] 

where C a (q±, q») is the normalized Fourier transform of the ath component of the magneti- 
zation density. In the simulations, we measured the time-evolution of \l/ and, after producing 
a rough estimate of their relaxation times, determined the time averages (&) and (\I/ 2 ) in the 
steady state. The runs typically ranged in length from 4 x 10 5 Monte Carlo sweeps (MCS) 
for the 12 x 9 systems to 4 x 10 6 MCS for the 32 x 64 systems. 

Figs. 2 and 3 display the results of our simulations. The data for \I/ (Fig. 2) clearly show 
a continuous transition to a phase-separated ordered state at /3± c ~ 0.68 (Note that the 
equilibrium Kosterlitz-Thouless transition |0| occurs at a lower temperature f3 c = 8/9 [0). 
The presence of a continuous transition is further supported by the cumulants for various 
finite-size systems (Fig. 3), g^ = 3(1 — (2/3)(^ 2 )/ (^))- In the limit of large sizes, they cross 



at the critical point j3^ c , as usual for continuous transitions [I5| . 

Since Figs. 2 and 3 show convincingly that the system does order, we now turn the 
question of how the two-temperature XY model can have an ordered state with long-range 
order. To address this issue we construct a field theoretic description of the model that 
indicates that the MW theorem fails because of the presence of effective long-range interac- 
tions generated by the two-temperature, diffusive dynamics. Although not rigorous, similar 
arguments have previously been used to discuss the generation of long-range interactions in 
Ising systems [|7|,|T6"1|, and in the spherical limit of the two-temperature diffusive dynamics 

We start with the 0(2) version of model B of Halprin and Hohenberg W71 which describes 



critical relaxation by Kawasaki dynamics towards the equilibrium of the coarse-grained XY 
model: 

d t $ = AV 2 [(-V 2 + t)$ + ^gU 2 } + V 

Here, 0(x, t) is the two-dimensional, coarse-grained order parameter, the parameters A and 
g are constants, and the relevant temperature dependence is contained in r, such that the 



critical point corresponds to r = 0. Furthermore, f/(x, t) is a Gaussian noise source which 
has zero mean and (r} a (x., t)r}p(x! , t')) = 2\5 a p5(t — t')V 2 5(x — x'). 

Now consider the generalization of model B to a rf-dimensional two-temperature model 
with spin exchanges that occur at one temperature in an m-dimensional 'parallel' subspace 
and at another temperature in the remaining d — m 'perpendicular' dimensions. In order 
to account for the different temperatures of the exchanges in the different subspaces of 
this model, the Laplacian operators and the noise term must be split into parallel and 
perpendicular parts. Furthermore, parameters such as r and g will have different values 
depending whether they are associated with the diffusion in the parallel or the perpendicular 
directions. For example, the model B term rV 2 will split into two terms: t±V 2 ± $ and T||<9 2 0, 
where d (Vj_) indicates a gradient over the m parallel (d — m perpendicular) directions. As 
r has now been split into parallel and perpendicular parts, the theory can describe various 
types of critical behavior depending on whether: (i) rj_ = and t\\ > 0, (ii) t\\ = and 
r_i_ > 0, or (iii) t± = t\\ = 0. The type of ordering ordering seen in Fig.l corresponds to (i) 
and in the following we restrict our consideration to that case. Keeping only those terms 
which are relevant (in the renormalization group sense) near the upper critical dimension 
d c = 4 — m, we arrive at the following equation 

d t $= AVi[(-Vi +r\\^ + r ± )j + ^gjtf] + V±- 

Note that only the perpendicular part of the noise, ff±, is relevant near d c . It has vanishing 
mean and (?7± a (x, t)r7_ L/3 (x / , t')) = 2\5 a p5(t - t')V±6(x - x'). 
Interestingly, this Langevin equation is of the form 

5(j) 

where fj satisfies the fluctuation dissipation theorem and therefore the above equation de- 
scribes the critical dynamics of an equilibrium system. The Hamiltonian, H, of that equi- 
librium system is most easily expressed in Fourier space as 

1 r k\ + T\\k\ + r ± k 2 , n I r - 

H= 2j k e ± -^M"*) + 4f J kl fc4 <P(kMk2)<P(ks)<P(h)5(ki + ■ ■ ■ + h), 
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which is the Hamiltonian of the the XY model with dipolar interactions. Thus, as in the 
Ising case [[7], we arrive at the conclusion that the mesoscopic critical properties of the 
two-temperature diffusive XY model are in the same universality class as an equilibrium 
model with long-range (dipole) interactions. Most importantly for the two-temperature XY 
model, however, we see the mechanism that apparently causes the Mermin- Wagner theorem 
to fail: the two-temperature diffusive dynamics generate effective long-range interactions in 
the steady state of the system. 

The static critical properties of the equilibrium XY model with dipole interactions have 
been studied using renormalization group methods for the m — 1 case |18| . The results show 
that the static structure function scales as 

S{k ± , k\\,r) = ^ 2+ri S{k ± /fi, fc||//i 1+A , r/fi 1 ^) 

which defines the exponents rj, v, and A. To second order in e = 3 — d, it has been found 
that v « 0.500 + O.lOOe + 0.054e 2 , rj « 0.0177e 2 , and A = 1 - rj/2. 

Fig. 4 shows the scaling plot of using the above values of the exponents with e — 1. 
The data collapse is not perfect, but it is consistent with the predictions of the field theory, 
considering that the field-theory predictions are only approximate, being to finite order in 
e. We caution, however, that a range of exponents [y = 0.6 ±0.1 and rj = 0.10 ± 0.15) also 
produce a collapse of the data. A more accurate determination of the value of the exponents 
would require data for a variety of system shapes and significantly larger system sizes than 
were used in the current study. 

In summary, our MC simulations taken together with the field-theoretic results suggests 
a simple picture: Long-range order does exist in the d = 2, two-temperature, diffusive 
XY model and it is caused by effective long-range interactions generated by the anisotropic 
diffusive dynamics coupled with the violation of the fluctuation-dissipation theorem. In light 
of this result, we return to the question posed at the outset: Can the LP argument or the 
MW theorem be extended to nonequilibrium steady states? In general, the answer to this 
question appears to be no, even for systems with purely local interactions and dynamics. 
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However, it still remains possible that the LP arguement or the MW theorem can be extended 
to nonequilibrium steady-state systems which evolve with some restricted class of dynamical 
rules. Unfortunatly, our knowledge in this field is rather limited; we know only that long- 
range order can appear in low- dimensional nonequilibrium systems whose dynamics is either 
long-ranged |J, or involves anisotropic diffusion jj. Clearly, much more work needs to be 
done to fully answer the question posed in this paper. 
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FIGURES 

FIG. 1. A typical ordered configuration with long-range order. The _L and || directions 
correspond to the directions of the _L and || bonds, respectively. The configuration shows a spin 
wave in the perpendicular direction. 

FIG. 2. Monte Carlo results for the order parameter showing a transition to long-range 
order at j3± c « 0.68. The size and shape (L± x Ly) of the systems is indicated in the legend. Error 
bars are much smaller than the symbol size. 

FIG. 3. Monte Carlo results for the cumulant gi defined in the text. The data for the different 
finite-size systems cross asymptotically at (3± c « 0.68, in agreement with the data for shown 
in Fig. 2. The size and shape, Lj_ x Lp of the systems is the same in Fig. 2. Error bars indicate 
one standard deviation of statistical errors. 

FIG. 4. Finite-size scaling of Monte Carlo results for the order parameter using 
(3± c = 0.68, v = 0.65, and r\ = 0.02. The size and shape, L± x Ly, of the systems is the same as in 
Fig. 2. 
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